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Abstract 

We investigate the properties of the amphtude induced by the anomaly. In a 
relatively high energy region those amplitudes are constructed by the vector meson 
poles and the anomaly terms, in which the anomaly terms can be essentially evaluated 
by the triangle quark graph. We pay our attention to the anomaly term and make 
intensive analysis of the existing experimental data, i.e., the electromagnetic vr'' and 
u) transition form factors. Our result shows that it is essential to use the constituent 
quark mass instead of the current quark mass in evaluating the anomaly term from 
the triangle graph. 



*Pcrmancnt address: Aichi University, Miyoshi Aichi 470-02, Japan 

'I'Fellow of the Japan Society for the Promotion of Science for Japanese Junior Scientists 



Much interests have been paid to the non-abehan anomaly because of its fundamental 
property of quantum gauge theories. These kinds of anomaly are known to be caused 
by the quantum effects of triangle fermion loops shown in Fig 0. Their amplitudes in 
the low energy limit is known to be independent either of the higher order corrections 
or of the internal fermion masses. On the other hand the asymptotic freedom of QCD 
facilitates us to predict the hadronic amplitudes at extremely high energy region, where 
again the triangle diagram becomes dominant. Indeed this triangle amplitude correctly 
reproduces the high energy behavior obtained from the operator product expansion (OPE) 
technique [0, |[, up to the ambiguity of the coefficient factor. 

The important observation here is that the amplitudes induced by the anomaly are 
constrained not only by the low energy theorem^, ^ ^] but also by the high energy behavior 
because of the asymptotically freedom, and therefore we may expect that the amplitudes 
induced by the anomaly are essentially controlled by the triangle quark graph. The most 
well-known example of the processes induced by the axial anomaly is vr*^ 77. The 
triangle quark amplitude reproduces the experimental data for this process excellently. 

If we want to apply the anomaly induced amplitudes to the relatively higher energy 
region, the non-anomalous contributions come in the process in addition to the anomaly 
terms. Then the QCD corrections to this triangle graph becomes important. Especially 
we have to take account of the following nonperturbative QCD effects: 1) chiral symmetry 
breaking (generation of constituent quark mass); 2) confinement effects; 3) hadron contri- 
butions. Our interest here is investigating the processes induced by the axial anomaly in 
the high energy region. 

In the separate paper 0, taking account of the above effects, we have proposed an in- 
terpolating formula for the amplitude induced by the axial anomaly. There the amplitudes 
are constructed by the "vector meson pole terms" and the "anomaly terms", the latter 
being essentially evaluated by the quark triangle graph of Fig. |I[ 

In this paper we pay our attention to the anomaly term and make intensive comparison 
with the existing experimental data. The candidates for investigating the structure of this 
triangle anomaly amplitudes are the 7r° transition form factors. 

The result shows that, in addition to the vector meson poles, it is essential to use the 
constituent quark mass instead of the current quark mass for the internal quark lines of 
the triangle graph. 
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Let us start with the following three-point function: 

TAb,,up^P^ q) = -^J d'xd'ye-^^ y^'^-^ (0| i:j^''{x)f'{Q)f,{y) |0) , (1) 

where q = k—p, is the axial- vector current which generally couples to Nambu-Goldstone 
(NG) boson P: = ipT^^lbi^, (0| 75(0) |P(k)) = ikjp, and j^^(x) and j^P{0) are the 
relevant currents, j^'^ = tpQ'^'y^tp with being the charge matrix of the quark field tp. 
Denote the amplitude the NG-boson pole removed as T^^''^'^P{p,q): 



I ■■ ■ [p, q) = I ■■ ■ [p, gj -h - 

n 

where 



TAb,,up^P^ g) = T^^'^^^(p, q) + ^i_^T^^.-P(p, q), (2) 

mi — k'^ 



T^^'^f>{p^ q) = -ij d'xe-''^^ (0| Tj^''ix)f^iO) |P(k)) . (3) 
If one fixes the current j^'^ = j^^ , the general form of this T^^'^^ij), q) is given by[0, 

X 



+ ip^q^q' - {p ■ q)q''p')H$\p\ g') + £"^>„g/JFp(p^ g^) 



(4) 



where A'c(= 3) is the number of colors and we use the electromagnetic current conservation: 
Pi,T^^'^P{p^q) = 0. We have normalized in such a way that Fp{p^ = 0, g^ = 0) = 1 (see 
Eg ■ (p!oD ) ■ If we further set = (real photon), we have the electromagnetic transition 
amplitude (P(k)| j^''(0) |7(p)), which is expressed in terms of the vector and axial form 
factors Pp(p^,g^) and F^\p'^,q^). In the case where the current j^''{0) couples to the 
photon or Z boson, the axial vector form factor vanishes, and only the term proportional 
to £°'t^'^p remains. Thus the only anomaly term Pp(p^, g^) contributes to such processes. In 
this paper, we consider this form factor Pp(p^, g^), which is directly related to the anomaly 
part of this three-point amplitude. 

To see how this form factor is incorporated to the anomaly term, we study the anoma- 
lous Ward-Takahashi identity expressed as[^ |TU| 

d^j^ = 2moj5 + -^FF^ (5) 

where js is the corresponding pseudoscalar density: 2moj5 = 2ijjM.T^i'j5i/j (A^: the mass 
matrix of the quarks). This leads us to 

k^T^B,,u,^p^ g) = M^^''^'^(p, q) + ^tT (t {q^, g^}) 5"^>„g^, (6) 
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where 



M^^^-P{p, q) = I d^a;Ae-*^"+*' (0| Tj^''{x)f^mmoMy) |0) • (7) 

This amphtude M^^'^^ip^ q), as well as the three-point function T^^'^'^'^^p, g), contains the 
NG-boson pole contributions. We extract the NG-boson pole contributions in both sides: 



fpk' 



rrip — k"^ 



+ 



4vr2 



tr(r{g^g^}). 



where use has been made of (0| 2moj5(0) |P) = rripfp. Thus the transition amplitude for 
NG-boson is obtained: 



r^B,up^p^ g) = - \k,f^'''^-'P{p, q) - M^'^'^'^ip, q) 
Jp '- 

which is rewritten in terms of the invariant amplitude 

N p2 

r^^^^P{p,q) 



tr(T{g^,Q^})£"^>,g^. 



(9) 



4vr2/p 



tr(T{Q^,Q^})£°^>„g^ 



l-np\q\k') 



(10) 



From this expression we easily see that T{p'^,q'^,k'^) corresponds to Fp{p'^,q'^) of Eq.(|^). 
It is well known that the triangle diagram containing fermion loops couples to vector or 
axial vector currents leads to anomalies, and in the standard renormalization procedure the 
gauge invariance guarantees that any higher order corrections do not modify the structure 
of this anomaly. 

So let us first calculate the quark triangle graph shown in Fig. |l| which is expected 
to give the dominant contribution to the function T(p^,g^,A;^) in Eq. (p!OD . The result is 
expressed asi 



np\q\e' 



2E^T,QfQf 



triangle tr (T {Q^, Q^}) 



[dz] 



mt 



with the Feynman parameter integral defined by J[dz] = 2 dzidz2dz^5{l — zi — Z2 — z^) 
It is easy to see that in the low energy limit in which k"^ = p^ = q^ = 0, this reduces to 



T(p2 = 0, g2 = 0, k'^ = 0) 



triangle 



(12) 



■l-Here we restrict ourselves to neutral currents. 
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Figure 1: The quark triangle anomaly graph 

independently of rrii. This is the expression that the anomaly term is exactly reproduced 
by the lowest order triangle graph. 

At this point, we should stress that the quark mass rrii appearing in T(p^, q'^,k'^) is 
taken to be the constituent quark masses instead of the original current quark masses. 
This reflects the most prominent feature of QCD: since the chiral symmetry has been 
spontaneously broken and there appear NG-bosons, the quarks necessarily acquire dynam- 
ical masses. This replacement never changes the value of the low energy limit, because 
T(p^ = 0, = 0, A;^ = 0) = 1 independently of the quark mass (see Eq.(|l2D). 

In the following we restrict ourselves to the ir^ case, where only the u and d quarks 
come into play. Then = Y^iqiTa^l^qu qi = {u,d), Ti = (1/2,-1/2), and 2moj5 = 
J2i 2?Tj,oj^jTji75gj. The triangle graph contribution to T(p^, g^, = 0) is obtained as 

/(p^gV) - Tip^q^e = 0) = / (13) 



where m (=m„ ~ m^) is the constituent quark mass (of the u and d quraks). 

If one photon is on its mass shell (e.g., p"^ = 0), then the function I{p'^,q'^;m) reduces 

to 



J{q^]m) = /(p^ = 0,g^;m 

,2 



m 

g2 



\/Am''^ — g2 — y/—q'^ 



(14) 



which gives the 7r°7 transition form factor. The high energy behavior of this function is 
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given by 

J{q ;m) ^ — 



(15) 



It is well known that, in the low energy limit, the higher order corrections do not 
change the value of the three-point function induced by the anomaly in Eq.(|I|), or T(p^ = 

= k"^ = 0) itself. However, QCD corrections do contribute appreciably to its higher 
energy behavior. Especially nonperturbative QCD effects are alleged to take place: 1) 
chiral symmetry breaking (generation of constituent quark mass); 2) confinement effects; 
3) hadron contributions. 

As we have mentioned below Eq.(|r^), we have already taken account of the effect 
1) by replacing the current quark masses by the constituent masses of the quarks. The 
next important QCD effects are represented as a rich hadron spectrum. As an average, 
most part of those resonances with broad widths may be already involved by the above 
replacement in the triangle contributions of Fig. 0. 

However, one should further modify the functions /(p^, g^; m) and J(g^; m) due to the 
effect 2). The graph in Fig. |l| shows the contribution of the qq threshold at = Am? 
[q^ = 4m^) from which an imaginary part emerges as shown in Fig. ^ (the dotted lines). 
However, since the quarks are confined, the intermediate states are not multi-quark states 
but actually multi-hadron states (27r, Air, ... etc.). The above effect may be properly taken 
into account by smearing out the original J(g^; m). Here we adopt the following function: 

~ 2 , (ln(l + x))' ln(l + x) , o^ n .ln(l + x) / A\ q^ 

J{q^;m) = — H ^ + ax exp{-f3x) + 27n^ ^ exp ] , x = 

X X X \ X J 

(16) 

Both J and J are shown in Fig. 0, where we take the parameter choice a = 1.4, (3 = 0.35 
and A = 2.5. From this figure, we see that the modified amplitudes J are properly smeared 
out. This form has been chosen so as to coincide with the original J (the dotted lines) 
both in the high and low energy limits. 

Finally, we must take account of the dominant hadron poles as the effect 3). Among 
hadron spectrum, NG-bosons, which have been put in as elementary fields, are of course 
most important. The next important poles are lowest energy bound states, i.e., the vector 
mesons. This can be done by introducing the vector meson poles to the amplitude, where 
we should be careful to include them to avoid the double counting and to satisfy the low 
energy theorem Eq.(^2|). 
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Real part of J{q^]m) and J{q^;m) 



1 1 


1 1 1 1 

Re [J] ■ ■ ■ - 
Re [J] ~ 






1 1 1 1 1 1 



10-2 10-^ 1 10^ 10^ 10^ 10"^ 10^ 



Imaginary part of J{q^;m) and J(q'^;m) 




10-2 10-1 1 IQi 102 IQS 10^ IQS 

q'^/m^ 



Figure 2: The QCD corrected-improved functions J{q^]m) (solid lines) comparing with 
the original functions J{(f'] m) (dotted lines), where we take the parameter choice a — 1.4, 
(5 = 0.35 and A = 2.5. 
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Taking all the above mentioned corrections into account, we propose a probable ex- 
pression of form factors. The 7r°7*7* transition form factor isi 

+ [[D,{p^) + D^{p')} J{q'- m) + {D,iq') + DUq')} Ap'; m) 
+ m [Dp{p')D^{q') + D,{q')D^{p^)] , (17) 

where Dp{p^) and Di^lp"^) are the Breit-Wigner type propagators for p and u mesons: 



with Fp and F^^ being the widths of p and u mesons, respectively. The weight of Dp and 
has been chosen to be consistent with the flavor symmetry. Since the above form 
factor should satisfy the low energy theorem, the normalization condition Eq.(|l2D gives 
Ki + 4^2 + /^s = 1 (note that DpijP' = 0) = D^^p"^ = 0) = 1). This expression is most 
elegantly derived in the framework of hidden local symmetry (HLS) (see, for a review, 
Ref . |TT| ) . In the following we use the parameters C3 and C4: C3 = 2/t3; C4 = 1 — ki, 



where C3 and C4 correspond to the terms appearing in the Lagrangian constructed in the 
framework of HLS. Then of course, this form factor again satisfies the low energy theorem: 
Ft,o^,^,{j)^ = 0, = 0) = 1, independently of the parameters C3 and C4. This together with 
the coefficient factor g-y^^, = — e^/(47r^/^), excellently reproduces the experimental value 
of the 7r° 77 process: this is nothing but the expression of the low energy theorem. 
In the high energy region, on the other hand, in addition to the J(g^;m) ~ (m^/g^) 
contribution, we have terms Dp{q^) ~ rn'j^/q'^ and D^^q"^) ~ m'^/q'^. If one notices that 
rrip ~ ~ 0{m), it is easy to see that p and uj pole contributions have the same high 
energy properties as I and J. 

Now, we compare the above proposed form factor with the experimental data and 
examine the q^ dependence of the anomaly terms I and J. To do this, we first determine 
the parameter C3 and C4 from the experiment. 

By setting = we have the vr" electromagnetic transition form factor: 

C3 + C4\ ~, 2 X C3 + C4 



F.o,{q') = 1 - ^ J{q';m) + ^ {Dp{q') + D^{q')} , (19) 



§The function / is obtained by the same procedure as J. 
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in which (cs + C4)/2 is to be determined. The comparison of our formula with the experi- 
mental data[|l^ indicates 

= 1.0. (20) 

This value implies that the vector meson dominance (VMD) is almost realized in this 
process. The result Eg . (pO]) implies that the dependence of J is not so important for 
the 7r°7 transition form factor, so far as we analyze the existing data (see Ref. 0). So we 
can not get the information of the function J(g^; m) from this process alone. 

On the other hand, the vector meson form factors can also be obtained from our 
formula. The expression for the cuvrfl transition form factor is given by extracting the 
cij-pole contributions from Eq.(^) and by making proper normalization of the amplitude 
we find 

(easel) F^{q^) = -cJ{q^-m) + {l + c)D,{q^), g = . (21) 

C3 + C4 

where c is to be determined. We take here m = mp/2 (case 1). For reference, we examine 
the case where we adopt the current quark mass m = mo — 5MeV in J(g^;mo) (case 2), 
and the case in which the triangle anomaly term is taken to be a constant (case 3): 

(case 2) Fi^-\q') = -dJ{q'; m,) + (1 + c)Dp{q') , (22) 

(case 3) F^{q^) = -g + (1 + c)Dp{q^) . (23) 

The form factor of case 3 was often used in the extensive analysis of the form factors in 
relatively lower q^ region (see, for example Ref. ||13||). 

Let us determine the value of c for each case using the experimental data r(u; — > 
TT^fi'^fi^). It is convenient to use the following expression: 



T{u ^ 7r07) J Ami 37r q^y^ q^ ) \\ q^ 

3/2 



q2 _ 



X 



q'^ \^ Am^a"^ 
1+ ^ ' 



m?, — mi j (m?, — mlV 



2 



FUqI , (24) 



where is the intermediate photon momentum (or invariant mass of final muons). The 
experimental data|]14| shows 



r(.^.VV-)^(9.6±2.3)xlO-^^ (25) 
r(c^ ^ TT^-f) (8.5 ± 0.5) X 10-2 ^ ' ' ^ ' 



^We could also obtain the p^Tr^ transition form factor. However, the branching ratio -B(/0° ^°7) 
(= (7.9 ± 2.0) X 10-**) is far smaller than B{lu —> tt'^^) (= (8.5 ± 0.5) x 10^^). Moreover, no experimental 
data for p° — > 7r°7* has been reported. 
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which give Jl 



(case 1) c = 0.49+[j-}| for F^{q^ 



(case 2) c = O.ll+jjj^, for F^™o)(g2), 

(case 3) c = 0.88+[}-{^, for F^(g^). (26) 

The average values for all the above cases show that the complete p meson dominance is 
incapable of describing the cuvr^ form factor (VMD corresponds to c = 0, see Eq.(pl|)). 

Now that we have found that the cjvr^ transition form factor receives appreciable contri- 
butions from the anomaly term, we use the uir^'j* process to get the information on the 
dependence of the anomaly term from the experimental data. The available experimental 
data for ujh^'j* lies in a relatively wide energy range (up to 1.4GeV). We use the value of 
c for each case obtained above, and compare the above three cases. The predicted curves 
are shown in Fig. \i (a) (b)| . In those figures the three data points around = 0.4GeV^ 
are located far from the theoretical curves. We should remark that the deviation of the 
above data points should not be taken seriously: because the above region corresponds to 
the kinematical boundary, the event rate is very rare. As a result, the behavior of the form 
factor strongly depends on the normalization uncertainty|T^]. From Fig. p (a) (b)] we find 
that the form factor F^{q^) (case 3) can not reproduce the high q^ data, the theoretical 
values are much larger than the data points. For the case 2 {F^°\q^)) the suppression of 
the form factor is too strong and the theoretical values are too small to fit the data points 
in the high q^ region. On the other hand, aside from the three points mentioned above, 
the improved form factor F^^^q^) (case 1) is found to be in good agreement with the data 
in the whole energy region. This successful result is owing to the triangle anomaly term 
described by the modified function J{q^\ m) in addition to the vector meson pole effect. In 
particular, in the above analysis we have taken the constituent quark mass m = mp/2. 

One might wonder whether the above results depend on the smearing procedure of 
the function J. To check this we show the curve obtained from the form factor with the 
unsmeared function J in Fig From this figure, we see that the smearing procedure does 
not affect the value of the form factor so far as in the region where the experimental data 
exists, although it does around the threshold (g^ ~ 0.6GeV^). 

In conclusion, the experimental data in a relatively wide energy range is found to be 
reproduced very well by using the QCD corrected triangle amplitude as the "anomaly 



"The another solutions for c, for example, c = —2.5 ±0.1 (case 1), are excluded by experiments. 
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term" . It must be stressed that the constituent quark mass plays an essential role in the 
triangle anomaly term. 

It is a common understanding that the low-energy effective Lagrangian for quarks and 
gluons is described by the chiral Lagrangian in terms of the NG-bosons. Their anomaly 
term has been introduced by the Wess-Zumino action[|l^, |1^ in the effective Lagrangian. 
Another way to incorporate the anomaly term is to introduce the quarks (to be more exact, 
we should say "constituent quarks") and gluon fields in place of the hadrons other than 
NG-bosons in the chiral Lagrangian [[T^ . 

At extremely low energy limit, the above two effective Lagrangians work very well for 
reproducing the phenomena and they are completely equivalent because the latter induces 
the anomaly term via loop effects. This is the result of low energy theorem. If one wants 
to apply the Lagrangian to nonzero but relatively lower energy scale, the vector mesons 
should be properly included in the Lagrangian. So far as we examine the processes in 
which the complete VMD is realized, the whole anomaly terms are replaced by the vector 
meson pole terms. Indeed we have seen that the anomaly term disappears from the 7r°7 
transition form factor and we cannot check the dependence of the anomaly term from 
this process. On the other hand, we have found that ujtc^ transition form factor, in which 
the complete VMD hypothesis is not valid, recieves appreciable contributions from the 
anomaly term. For such processes the above two approaches give different predictions. 
Our result indicates that the WZ action alone cannot reproduce the experimental data 
for the uj7T^ transition form factor in the relatively high enenrgy region, because the WZ 
action gives the constant anomaly term. This is in remarkable contrast to the anomaly 
term calculated from the QCD corrected quark triangle graph. In this sense our treatment 
here may have somewhat similar spirit as the latter approach. 

Also the high energy behavior of the anomaly term, with the usage of the constituent 
quarks in the internal loop, seems to suggest the possibility that this anomaly term can be 
applicable even to the phenomena in the extremely high energy region, say for examples, 
Z — > P7 (g^ = m|) or W decay processes]^, |^, ||, ^ |T^. The processes induced by the 
axial anomaly may yet provide us with the information on how to construct the effective 
theory applicable to hegher energy regions. 

We would like to thank Taichiro Kugo for comments and instructive discussions about 
the property of the anomaly. 
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Figure 3 (a) (b): The form factors in the energy range below 1.4 GeV: ( 1) the 

improved form factor F^^q"^); (b: case energy data (Exp.l) is given 

in Ref.|15|, and the high energy data (Exp. 2) is translated from the cross section data|jl9 
The values of c are determined using T{uj vr^/i+zi") (see Eq.(|23)). 
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(c: case 3) 
"1 1 r 



1 1 1 

c= 1.04 — 
c = 0.88 — 
c = 0.70 ■ ■ ■ - 

Exp.l 

Exp.2 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

g2 [GeV2] 

Figure 3 (c): The form factors in the energy range below 1.4 GeV: (c: case 3) F^{q^). 



2\|2 



100 




q' [GeV^] 



Figure 4: The form factor in the energy range below 1.4 GeV using the unsmeared function 
J{q^]m). 
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